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in ■ Abstract 

We elaborate on the proposal that the observed acceleration of the Universe is the result of 

in . 

the backreaction of cosmological perturbations, rather than the effect of a negative-pressure dark- 
energy fluid or a modification of general relativity. Through the effective Friedmann equations 



describing an inhomogeneous Universe after smoothing, we demonstrate that acceleration in our 
local Hubble patch is possible even if fluid elements do not individually undergo accelerated ex- 
pansion. This invalidates the no-go theorem that there can be no acceleration in our local Hubble 
patch if the Universe only contains irrotational dust. We then study perturbatively the time behav- 
ior of general-relativistic cosmological perturbations, applying, where possible, the renormalization 
group to regularize the dynamics. We show that an instability occurs in the perturbative expansion 
involving sub-Hubble modes. Whether this is an indication that acceleration in our Hubble patch 
originates from the backreaction of cosmological perturbations on observable scales requires a fully 
non-perturbative approach. 
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I. INTRODUCTION 



Recent observations of the expansion history of the Universe indicate that the Universe 
is presently undergoing a phase of accelerated expansion . The accelerated expansion 
is usually interpreted as evidence either for a "dark energy" (DE) component to the mass- 
energy density of the Universe, or for a modification of gravity at large distances. In this 
paper we explore another possibility, namely that the accelerated expansion is due to the 
presence of inhomogeneities in the Universe. 

In the homogeneous, isotropic, Friedmann- Robert son- Walker (FRW) cosmology, the ac- 
celeration (or deceleration) of the expansion may be expressed in terms of a dimensionless 
parameter q, proportional to the second time derivative of the cosmic scale factor a. It 
is uniquely determined in terms of the relative densities and the equations of state of the 
various fluids by (overdots denote time derivatives), 

^ _ _aa _ q total + ^ Wi rii, (1) 
a i i i 

where ^total is the total density parameter and the factors fli are the relative contributions 
of the various components of the energy density with equation of state Wi = Pi/ Pi (Pi and 
Pi being the pressure and energy density of 2-th fluid). The expansion accelerates if q < 0. 
Observations seem to require DE with present values w^e ~ — 1 an d VL DE ~ 0.7 Q|. The 
negative value of w^e, indicating a violation of the energy condition w > —1/3 j^j], is 
usually interpreted as the effect of a mysterious dark energy fluid of unknown nature or a 
cosmological constant of surprisingly small magnitude. 

The existence of a negative-pressure fluid or a cosmological constant would have profound 
implications for physics as well as cosmology. While the observational evidence for the 
acceleration of the Universe is now compelling, it is important to keep in mind that the 
evidence for dark energy is indirect; it is inferred from the observed time evolution of the 
expansion rate of the Universe. What is known is that the expansion history of the Universe 
is not described by the expansion history of an Einstein-de Sitter Universe (a spatially flat, 
matter-dominated FRW model). While such a departure may be caused by dark energy, 
there are other possibilities. One possibility is that general relativity is not a good description 
of gravity on large distance scales. Another possibility is that the Universe is matter- 
dominated and described by general relativity, and the departure of the expansion rate from 
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the Einstein-de Sitter model is the result of back reactions of cosmological perturbations. 
This explanation is the most conservative, since it assumes neither a cosmological constant, 
a negative-pressure fluid, nor a modification of general relativity. 

In this paper we explore the possibility that backreactions of cosmological perturbations 



is the source of the accelerated expansion 



12. 131. The idea is as follows. 



We know there exist cosmological perturbations; after all, the Universe is inhomogeneous. 
To describe the time evolution of a patch of the Universe as large as our local Hubble radius 
one has to construct the effective dynamics from which observable average properties can be 
inferred. Of course, this implies a scale-dependent description of inhomogeneities. Suppose 
further that our Universe is filled with pressureless matter and no DE. If inhomogeneities 
evolve with time, a local observer would infer that our Universe is not expanding as a 
homogeneous and isotropic FRW matter-dominated Universe with Hubble rate Hit) oc t 2//3 , 
where t is cosmic time. On the contrary, the Universe would appear to have an expansion rate 
with a time evolution that depends on the nature and time evolution of these perturbations. 
Potentially, this could lead to an accelerated expansion. 

Our idea is actually intimately connected with the general problem of how the nonlinear 
dynamics of cosmological perturbations on small scales may affect the large-scale "back- 
ground" geometry. Let us start by discussing this issue in some generality. 

The standard approach to cosmology is based both on observational facts, such as the 
near-perfect isotropy of the Cosmic Microwave Background (CMB) radiation, and on an a 
priori philosophical assumption, the so-called Cosmological or Copernican Principle. Ac- 
cording to the Cosmological Principle all comoving cosmic observers at a given cosmic time 
should see identical properties around them (isotropy around all cosmic observers implies 
homogeneity, hence the FRW line element). The Cosmological Principle allows one to cir- 
cumvent our inability to obtain information about the Universe outside our past light-cone 
by assuming that a symmetry principle is valid everywhere. By using the Cosmological 
Principle, we assume that we are able to determine conditions many Hubble radii away from 
us by using observational data within o ur p ast light-cone, whose region of influence is, by 
definition, limited to one Hubble radius [14]. 

An alternative procedure, dubbed Observational Cosmology, has also been proposed. 
It aims at constructing a cosmological model solely in terms of observational facts, thereby 
avoiding any a priori assumptions of global symmetry. It dates back to the works by Kristian 
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and Sachs in 1966 and Ellis in 1984 |"L4j ]. A remarkable feature of this approach is that, 
by using Einstein's equations, the dataset observable within our past light-cone is precisely 
sufficient to determine the space-time and its matter content within the same light-cone (see 
Ref. Q and references therein). 

A crucial ingredient of the Observational Cosmology approach, which is shared by any 
realistic cosmological model-fitting procedure, is smoothing. Observations tell us that the 
Universe is far from homogeneous and isotropic on small scales. Observationally, we know 
that homogeneity, e.g., in the galaxy distribution, is only achieved over some large smoothing 
scale (see e.g., Ref. (l6j). When we refer to homogeneity and isotropy of the Universe we 
tacitly assume that spatial smoothing over some suitably large filtering scale has been applied 
so that fine-grained details can be ignored (see in this respect the discussion in Refs. }l7, rsl]). 
In other words, by the mere assumption that the same background model can be used to 
describe the properties of nearby and very distant objects in the Universe, the smoothing 
process is implicit in the way we fit a FRW model to observations. Cosmological parameters 
like the Hubble expansion rate or the energy density of the various cosmic components are to 
)e considered as volume averaged quantities: only these can be compared with observations 

There is, however, a technical difficulty inherent in any smoothing procedure. While mat- 
ter smoothing is somewhat straightforward (e.g., in the fluid description), smoothing of the 
space-time metric is more complex and immediately leads to an important and unexpected 
feature, pointed out by Ellis [IJ]. Let us assume that Einstein's equations hold on some 
suitably small scale where the Universe is highly inhomogeneous and anisotropic. Next, 
suppose we smooth over some larger scale. Einstein's equations are nonlinear: smoothing 
and evolution (i.e., going to the field equations) will not commute. Hence, the Einstein 
tensor computed from the smoothed metric would generally differ from that computed from 
the smoothed stress-energy tensor. The difference is a tensor appearing on the right-hand 
side (RHS) of Einstein's equations that leads to an extra term in the effective Friedmann 
equations describing the dynamics of the smoothing domain. 

How can this fact be related to the acceleration of the Universe? The answer is that 
this extra source term need not satisfy the usual energy conditions (according to which our 
Universe can only decelerate) even if the original matter stress-energy tensor does. The fact 
that the effective stress-energy tensor emerging after smoothing could lead to a violation 
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of the energy conditions was originally recognized by Ellis Q] 1 As we will discuss in 
Sec. HH explicit calculations of the effective Friedmann equations [3, 



21 



22 confirm that 



acceleration is indeed possible even if our Universe is filled solely with matter. 

A closely related question is what is the appropriate scale for which the smoothing pro- 
cedure can fit the standard picture of a homogeneous and isotropic Universe on large scales? 
Our choice will be that of smoothing over a volume of size comparable with present-day 
Hubble volume. The precise size of the averaging volume does not matter, provided it is 
large enough that the fair sample hypothesis applies, i.e., that volume averages yield an 
accurate approximation of statistical ensemble averages. We will nonetheless refer to scales 
within (outside) the averaging domain as "sub- Hubble" ( "super-Hubble" ) perturbations. 
In doing this we are however promoting our super-Hubble, or "zero" -mode, to the role of 
FRW-like background. 

The next question is what are the scales that determine the dynamics of our local back- 
ground. To answer this question we have to recall what happens in the standard FRW 
models. The evolution of the global scale factor a(t), the zero-mode of FRW models, is 
fully determined by the matter content of the Universe through the value of ^total an d 
Qj, and via the equation of state to; of its components. That is where microphysics enters 
the game. In other words, in the standard FRW picture the evolution of the Universe as 
a whole is determined by the dynamics of matter on sub-Hubble scales. Similarly, in the 
Observational Cosmology approach the dynamics of our local background must be deter- 
mined by the observed behavior of matter inhomogeneities within our past-light-cone. That 
is where the backreaction of sub-Hubble inhomogeneities enters the game. This picture will 
become clear in Sec. |Hj where we will introduce two scalars, the so-called kinematical back- 
reaction Qd and the mean spatial curvature (R)d, that enter the expression for the energy 
density and pressure in the effective Friedmann equations governing the mean evolution of 
our local domain D. The crucial point is that in the fully general relativistic framework 
these two scalars are linked together by an integrability condition (which has no analogue 
in the Newtonian context), whose solution provides the effective equation of state of the 
backreaction. In order to solve this equation and establish the typical size of these terms 
one needs a non-perturbative and non-Newtonian approach to the evolution of cosmological 



1 A closely related discussion can also be found in Ref. [2o| . 
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irregularities, as pointed out in Refs. 



23 



24|. 



The fact that the average dynamics naturally leads to new terms in the source implies that 
it is legitimate to use an effective Friedmann description, provided one takes into account 
that the effective sources of these equations contain the back-reaction terms. 

What will result from our analysis is that the evolution of sub-Hubble perturbations leads 
to an instability of the perturbative expansion due to the presence of large contributions 
which depend on a combination of Newtonian and post-Newtonian terms. This instability 
indicates that the effective scale factor describing the dynamics of our local Hubble patch is 
fed by the evolution of inhomogeneities within the Hubble radius. This cross-talk between 
the small-scale dynamics and the effective average dynamics described by super-Hubble, 
or "zero"-mode, playing the role of FRW-like background might be the crucial ingredient 
of backreaction that can lead to cosmic acceleration without dark energy. In Ref. a 
deviation from the pure matter-dominated evolution was obtained by a combination of sub- 
and super-Hubble modes generated by inflation, the latter being improperly used to amplify 
the backreaction. In this paper, we will show that the deviation from a matter-dominated 
background is entirely due to the nonlinear evolution of sub-Hubble modes which may cause 
a large backreaction (technically due to the disappearance of the filter modeling the volume 
average), while the super-Hubble modes play no dynamical role. 

At this point it may be useful to contrast the differences between our approach dealing 
with inhomogeneities with the traditional approach. In the traditional approach one averages 
over inhomogeneities and forms a time-dependent average energy density (p(x, t)) (although 
the standard procedure is to calculate averages with the unperturbed spatial metric!). One 
then uses for the dynamics of the "zero- mode" [a(t)] the dynamics of a homogeneous universe 
with energy density p(t) = (p(x,t)). One then regards inhomogeneities as a purely "local" 
effect, for instance, leading to peculiar velocities. In this approach inhomogeneities can not 
result in acceleration. 

In our approach, we take note of the fact that the expansion rate of an inhomogeneous 
universe of average density {p(x, t)) (using the inhomogeneous spatial metric to calculate the 
spatial average!) is not the same as the expansion rate of a universe with the same average 
density. In order to account for this we encode the expansion dynamics into a new zero 
mode (or scale factor) a^(t) (which will be properly defined in the next section). It is the 
dynamics of this renormalized scale factor that is best used to calculate observables and will 
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determine whether the Universe accelerates. In our approach the effect of short-wavelength 
inhomogeneities is not just a local effect, but renormalizes the long-wavelength dynamics. 

Our paper is organized as follows. In Sec. |H] we summarize the effective Friedmann 
description of an inhomogeneous Universe after smoothing. In Sec. Illll we discuss how accel- 
eration in our Hubble patch can result from the backreaction of perturbations. Conclusions 
are drawn in Sec. IIV1 The Appendix presents the main results of a fourth-order gradient- 
expansion technique. 



II. EFFECTIVE FRIEDMANN EQUATIONS IN AN INHOMOGENEOUS UNI- 
VERSE 

The goal of this section is to compute the time dependence of the local expansion rate 
of the Universe. For a generic fluid we may take the four- velocity to be = (1,0), which 
amounts to saying that a local observer is comoving with the energy flow of the fluid. For 
the case of irrotational dust considered in this paper we have the freedom to work in the 
synchronous and comoving gauge with line element 

ds 2 = -dt 2 + hij(x, t)dx l dx j , (2) 

where t is cosmic time. 

A fundamental quantity in our analysis is the velocity gradient tensor, which is purely 
spatial and symmetric because of irrotationality. It is defined as 

"', ^i>i, r (3) 

Here dots denote derivatives with respect to cosmic time. The tensor 0* ■, represents the 
extrinsic curvature of the spatial hypersurfaces orthogonal to the fluid flow. It may be 
written as 

Q i j = QS i j + <T i j . (4) 

Here is called the volume- expansion scalar, reducing to 3H (H is the usual Hubble rate) 
in the homogeneous and isotropic FRW case. The traceless tensor er* . is called the shear. 

The evolution equations for the expansion and the shear come from the space-space 
components of Einstein's equations (see e.g., Ref. |2i 



]). They read, respectively, (p is the 
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mass density, R and R l j are the spatial Ricci scalar and tensor respectively of comoving 
space-like hypersurfaces) 



& + & 2 + R = 12vrGp, (5) 
1 
3 



+ 6aV + R } - -R5^ = 0. (6) 



The — component of Einstein's equations is the energy constraint 

-Q 2 -2a 2 + R= 16nGp, (7) 
3 

where a 2 = \o % j<j 3 { . The — % components yield the momentum constraint 

ai j\i - |ej = 0, (8) 

where the vertical bar denotes covariant differentiation in the three-space with metric hij. 
The mass density, in turn, can be obtained from the continuity equation 

p = -ep, (9) 

whose solution reads 

P = Po (h/h y 1/2 , (10) 

where h = det hy and the initial conditions have been arbitrarily set at the present time 
to- Combining the expansion evolution equation with the energy constraint gives the Ray- 
chaudhuri equation, 

Q + ^Q 2 + 2a 2 + 4nGp = 0. (11) 

From the latter equation it is straightforward to verify that irrotational pressure fluid ele- 
ments cannot locally undergo accelerated expansion. (This point was emphasized by Hirata 
and Seljak Indeed, defining a local deceleration parameter and using the Raychaudhuri 
equation, one finds 

q = - (36 + 6 2 ) /0 2 = 6(a 2 + 2nGp)/Q 2 > 0. (12) 



While it is true that locally the expansion does not accelerate, it is incorrect to assume 
that acceleration can not occur when the fluid is coarse-grained over a finite domain. The 
reason is trivial: the time derivative of the average of G and the average of the time derivative 
of are not the same because of the time dependence of the coarse-graining volume. 



Let us denote the coarse-grained value of a quantity T by its average over a spatial 
domain D: 2 

We will take the domain to be comparable with the size of the present Hubble volume 3 . 
A first important property follows directly from the smoothing procedure itself: for a 




generic function T one has |21l . 12^ 

{F)n ~ (f)D = (F&)d ~ {®)d{T)d, (14) 

where we have not considered terms originating from the peculiar motion of the boundary, 
since we will eventually consider only comoving domains in what follows. In particular, for 
the local expansion rate one finds 

(&)d = (&)d + (& 2 }d-(&} 2 d > (Q)d- (15) 

Although (Q)d < — |(@ 2 )d < 0, the coarse-grained deceleration parameter qr> = 
—3(0)^/(0)1, — 1 is related to (Q) D , which is not the same as (0)d- It is precisely this 
commutation rule that allows for the possibility of acceleration in our local patch in spite of 
the fact that fluid elements cannot individually undergo accelerated expansion. This simple 

nn n 

argument circumvents the no-go theorem adopted in Refs. |26j, |2j| (and later in |28l|). ac- 
cording to which there can be no acceleration in our local Hubble patch if the Universe only 
contains irrotational dust. 

Indeed, let us follow the work of Buchert (21L 1^ and define a dimensionless scale factor 

fy \l/3 

a D (t)= -° ; V D = Vhd 3 x, (16) 



y D j ' Jd 

where Vp is the volume of our coarse-graining domain (the subscript "0" denotes the present 
time). As an averaging volume we may take a large comoving domain, so that D is constant 



2 Notice that one is not allowed to define the mean cosmological parameters only through an average 
over directions [3, [2^ as cosmological observables, such as the Hubble rate, receive unacceptably large 
corrections from the same Newtonian terms which become harmless surface terms when averaging over a 
large volume Q. We acknowledge discussions with U. Seljak about this issue. 
The correct definition of our comoving Hubble radius is R H (t ) = e~* M / die*'®. 
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in time. Alternative choices are however possible (see, e.g., Ref. |29| for a discussion of 
different averaging procedures). 4 

The coarse-grained Hubble rate H D will be 

H D = ^= 1 -(Q) D . (17) 

Old 6 

Notice that with such a coarse-graining, Hp coincides with the quantity H defined in Ref. 
^1). By properly smoothing Einstein's equations over the volume Vd, one obtains 2ll. I22I] 

^ = -^(Pc ff + 3P eff ), (18) 
fd D \ 2 8nG 

UJ = — p * (19) 

where we have defined effective energy density and pressure terms 

PcS = {P)D -16^-16^ (20) 

P _ Qd I ( R )d /on 
Pcff " " 16^G + (21) 

and we have introduced the kinematical backreaction 

Q D = 2 -((e 2 ) D -(&} D )-2(a 2 } D . (22) 

From the effective Friedmann equations, Eqs. (|THj) and fj!9jl . obtained by Buchert in Ref. 
(2]], one immediately obtains the continuity equation for our effective fluid 

p eff = -3H D (p cff + P cff ) . (23) 

Note that the smoothed continuity equation differs from the local continuity equation. On 
the other hand, owing to the fact that our coarse-graining volume is comoving with the mass 
flow, mass conservation is preserved by the smoothing procedure, implying 

(p) D = -3H D (p) D . (24) 

The two quantities Qd and (R)d are not independent. This can be seen by taking the 
derivative of Eq. (fTU|) and using Eq. The consistency of the system of Eqs. (fTHjl. (fT^j) . 

and (}24"|) requires that Qd and (R)d satisfy the integrability condition 

(a%Q D )' + a D (a 2 D (R) D )' = 0. (25) 



4 One could alternatively average over a volume of size of the order of the instantaneous particle horizon. 
The effective dynamics in this case will be accompanied by a stochastic source originated by the statistical 
nature of the perturbations. 
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One should stress that the latter equation, i.e., the link between kinematical backreac- 
tion Qd and mean curvature (R)d, is a genuine General Relativistic (GR) effect, having no 
analogue in Newtonian theory, as the curvature R of comoving hypersurfaces vanishes iden- 

implying that there exist globally flat Eulerian 



tically in the Newtonian limit 
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coordinates X 1 . Indeed, in the Newtonian case, it is immediate to verify that is exactly 
(i.e., at any order in perturbation theory) given by the volume integral of a total-derivative 
term in Eulerian coordinates 2^|, 

gNcwtonian = < V ■ [ll (V • ll) - (ll • V) u]) D , (26) 

where u is the peculiar velocity, so that by the Gauss theorem it can be transformed into a 
pure boundary term. It is precisely by this reason that any analysis of backreaction based 



31], 



is not 



on the Newtonian approximation, such as that recently performed in vl of Ref. 
relevant: it will invariably lead to a tiny effect, and to the absence of any acceleration. 
Indeed, if inhomogeneities only exist on scales much smaller than our Hubble radius and 
if peculiar velocities are small on the boundary of our Hubble patch then, within the New- 
tonian approximation, the standard FRW matter-dominated model can be applied without 
any substantial correction from the backreaction 1231. Such a drawback of the Newtonian 

fin 

approximation was also noticed in Refs. |6j, [7[ . This exact result demonstrates that in order 
to deal with the backreaction, going beyond the Newtonian approximation is mandatory, 
as also stressed in Ref. Studies of the average dynamics including the lowest post- 

Newtonian gradient terms in the weak field limit were considered in Refs. (3^ and v2 of 
3Jj. However, further and more sizeable terms are expected to contribute to the backreac- 
tion once the effective dynamics of the system (including the kinematical backreaction) is 
considered. We will come back to this issue in subsection IIIC. 

The GR integrability condition makes it clear how acceleration in our local Hubble patch 
is possible. Indeed, it is immediate to realize that the general condition for acceleration in 
a domain with mean density (p) D is 

Q D >AixG{p) D . (27) 
Moreover, a particular solution of the integrability condition for constant Qd and (R)d is 

Q D = ~(R) D = const., (28) 
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which, for negative mean curvature mimics a cosmological constant, A c g = Qd- More in 
general, if is positive, it may mimic a dynamical dark energy or quintessence. 

So far the considerations have been rather general. Now we write the spatial metric in 
the general form 



where a(t) oc t 2 l z is the usual FRW scale-factor for a flat, matter-dominated Universe and 
the traceless tensor Xij contains the remaining modes of the metric, namely one more scalar, 
as well as vector and tensor modes. 5 Next, when we consider the expansion in some domain 
D, we can split the gravitational potential \1/ into two parts: ^ = + where ^£ is the 
long-wavelength mode and ty s is a collection of short-wavelength modes. Of course "long" 
and "short" describe wavelengths compared to the size of the domain D. We can easily take 
into account the effect of ^ by noting that within D the factor e" 2 ^ 1 is a space-independent 
conformal rescaling of the spatial metric. Let us then write 



h .. = a 2(£) e -2*(x,t) [ < J i . +Xi .( X , i )], 



(29) 



hij = a 2 (t)e 



^ ^ hij (x, t) , 



(30) 



with = e 2 * s ( x '*) [§ij + Xij(x, t)]. The expansion scalar and shear then become 



6 



3_ + -3^, 



a 




(31) 



where 6 and a 1 ■ are calculated with hij. Note that and a 1 ■ do not depend explicitly on 
It should be kept in mind that the local expansion rate is 6, not 6. 
The Ricci tensor of comoving space-like hypersurfaces is given by 



K 3 = a 



(32) 



where i2 l • is the Ricci scalar of the metric h^ and the symbol V, denotes the covariant 
derivative in the 3-space with metric h^. For the Ricci scalar we find 



R = a - 2 e 2 * e [R + 4V 2 ^ - 2V fc tf<V fc tf *] (33) 
(R) D = a~ 2 e 2 ^ (R + 4V 2 ^, - 2V fc ^V fc ^) D . (34) 



5 Indices of Xij be raised by the Kronecker symbol: x l j = $ tk Xkj and x y = <? fe <P' 
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Notice that cld coincides with the scale factor a adopted in Ref. jsj], provided we take 

a D (t) =o(<)e-* l(t)+ *» (35) 

with 

^(Xobs, t) = In a — -In ^ -\//i(i 3 a;^ + const., (36) 

where the residual dependence of on the spatial coordinate x obs labels the individual 
comoving volume patch, i.e., the specific cosmic observer we are considering. Using Eq. 
(fTUj) we can rewrite Eq. in the form 

* e {t) = --ln((l + W) -1 ) An + const., (37) 

where <5frw is the density contrast with respect to the mean density of a flat, matter- 
dominated FRW (Einstein-de Sitter) model, defined through p = (1 + <5frw) / (G^Gt 2 ), and 
"in" denotes the initial time, which for simplicity we have taken to coincide with the end of 
inflation. For any quantity T 

(J~)d- ~ ~ — ^ X (38) 
Id \fh\rid^x 

By inspecting Eq. (|3*7j) one immediately realizes that acceleration may be achieved in those 
Hubble patches where the mean rarefaction factor ((1 + £frw) _1 ) An g r °ws fast enough to 
compensate for the Einstein-de Sitter expansion rate H = 2/3t. Note also that the integral 
defining tyg is dominated by the dynamics of the most underdense fluid elements, not by the 
densest ones, so the complex dynamics of highly nonlinear mass concentrations never enters 
the calculation; by the same reasoning, any intrinsic limitation related to caustic formation 
would not affect the validity of our backreaction treatment. 

The kinematical backreaction Qd is non- vanishing and gets contributions only from G 
and from the shear a 1 .•: 

Qd = 1(Q 2 )d-2(v 2 )d, (39) 

where we used the fact that (0)d = by construction. 

In order to have a qualitative understanding of why acceleration can be the natural 
outcome of the backreaction, let us rewrite the mean expansion rate in terms of the peculiar 
expansion rate 9, defined by 6 = 3H + 9: 

H - 2 I 1 ((* + <Wr X A) A» m 
D 3t 3 ((i + (W )-i) An ' 1 > 
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which shows that the mean expansion rate receives a correction with respect to the FRW 
background value by the peculiar expansion rate of mostly underdense regions (where 9 is 
largest), which give the largest contribution to the average. However, as we already noticed, 
acceleration may be achieved when Qd is positive and large enough. This requires a large 
variance of the volume expansion rate within the averaging domain. At early times, when 
perturbations are linear, G is narrowly peaked around its mean background value 3H. When 
non-linearities set in, the variance increases because of the simultaneous presence of largely 
under- and over-dense regions (in fact, counting only under-dense regions would reduce the 
variance leading to an under-estimate of Qd) [33]. 



III. THE APPEARANCE OF INSTABILITIES 



In order to discuss the dynamics of our local Hubble patch, one may proceed along two 
complementary directions. Either one tries to encode the effect of perturbations into the local 
scale factor ao as done in Ref. p, or one may try to construct an effective equation of state 
by computing the backreaction terms Qd and (R)r>. We will try to follow a combination of 
them to see under which circumstances acceleration in our Hubble patch can be achieved. 



A. The effect of super-Hubble modes 

The mean curvature generally depends on both sub- and super-Hubble modes. Since Qd 
does not depend explicitly on perturbations with wavelengths larger than the Hubble radius, 
one immediately concludes that if one considers super-Hubble modes only, Q D vanishes and 
from the integrability condition (R)d scales like a D 2 . Therefore the effect of pure super- 
Hubble perturbations is limited to generating a true local curvature term which may be 
important but can not accelerate the expansion of the Universe. 

As we already anticipated in the Introduction, the effective scale factor describing the 
dynamics of our local Hubble patch is fed by the evolution of inhomogeneities within the 
Hubble radius. This cross-talk between the small-scale dynamics and the effective average 
dynamics is a crucial ingredient, as also pointed out in Ref. 0. Before coming to this 
crucial point, let us pause for a moment and show that there is another more technical way 
to achieve the same conclusion about the role played by super-Hubble modes starting from 
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the spatial-gradient expansion of Einstein equations. 

The spatial-gradient expansion is a nonlinear approximation method suitable to describe 
the long-wavelength part of inhomogeneities in the Universe. This scheme is based on the 
assumption that observables like the local curvature can be expanded in powers of gradients 
of the perturbations. To account for the effect of super-Hubble modes at late times we 
may adopt the so-called renormalization group method applied to the gradient expansion 
of Einstein equations 2^1 ■ This will result in the renormalized long- wavelength solution to 
Einstein equations, valid also at late times until the long-wavelength perturbations enter the 
horizon. 

Here we sketch a non-perturbative technique to solve Einstein's equations in an inhomo- 
geneous Universe. A more detailed presentation of the method will be presented elsewhere 

Our approach makes use of a gradient-expansion approximation (see Refs. \3a\3A \M 

, |4Q, |41|, |42|). The idea is to describe the dynamics of irregularities in a Universe which 
contains inhomogeneities on scales larger than the Hubble radius. Working in the syn- 
chronous gauge one expands Einstein's equations starting from a space-dependent "seed" 
metric. The lowest order solution corresponds to the so-called long-wavelength approxima- 
tion, while adding higher-order gradients leads to a more accurate solution, which hopefully 
converges toward the exact one. 

The gradient-expansion technique amounts to keeping a finite number of spatial deriva- 
tives. This approximation technique is non-perturbative in the sense that by solving for 
the metric coefficients \I> and Xij U P to 2n spatial gradients one obtains terms of any order 
in the conventional perturbative expansion containing up to 2n gradients. For the purpose 
of this paper, working to four derivatives in \1/ and Xij w iU suffice. Note that because the 
scalar \1/ appears in the argument of an exponential in the way we write the spatial metric, 
our gradient-expansion sensibly differs from that used in Refs. 0, 0. 
Even if \1/ is obtained up to a finite number of spatial gradients, will necessarily contain 
gradient terms of any order. 

Since the cosmological perturbations are generated during inflation, it is natural to set 
initial conditions for the gravitational perturbations \I/ and Xij & t t ne en d of inflation (ef- 



fectively coinciding with t = t in = 0). If so, the spatial metric in t 
is given by hij = a 2 e~ 2 ^5ij where ( is the curvature perturbation 



re super-horizon regime 
It is related to the 
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so-called peculiar gravitational potential y(x) defined by ip = 3£/5, in a matter-dominated 
Universe. The comoving curvature perturbation is constant on super-horizon scales (up to 
gradients) when only adiabatic modes are present and the decaying mode is disregarded. 
Therefore, the initial conditions at t — are ^/j n = ^f(t = 0) = 5<p/3 and Xij(t = 0) = 0. 
Notice that since cosmological perturbations generated during single-field models of inflation 
are very nearly Gaussian with a nearly flat spectrum j^j], we infer that ip should be regarded 
as a nearly scale-invariant, quasi- Gaussian random field. 

From this analysis we see that \I/ contains at least a zero- derivative term: this will be 
our seed metric perturbation, which is the necessary ingredient for gravitational instability 
to develop. The traceless tensor Xij has at least two spatial gradients. The only exception 
might come from linear or higher-order tensor modes appearing at the end of inflation. 
Nonetheless, accounting for these contributions does not quantitatively affect our results. 

Let us adopt the gradient-expansion in order to obtain \1/ and xl including up to four 
gradients of the initial seed potential <p. Under this assumption, the inverse spatial metric 
can be written as 

h ij = a~V* (** - X ij + X\x kj ) ■ (41) 
Writing = 3H + 9 where H = a/a, one finds (up to higher-derivative terms) 

= -3i - X -x kl Xki (42) 

= \x l , - \x ik Xki + \x lk XmVy (43) 

To solve for these quantities one also needs the 3D Ricci tensor and Ricci scalar of the 
constant-time hypersurfaces. We start by calculating the 3D Christoffel symbols, which read 
(up to higher- derivative terms) 

T) k = , ^ - 8\ + 5 jk + \ (x l hk + X\j ~ XjS) + *'* Xjh ~ X%k- (44) 
The Ricci tensor and Ricci scalar read, respectively, 



R' 



K\ e 



3 a? a? 



2* 

+ V 2 *5* • + - (W) 2 <f ■ - x ik *,kj ~ x ik ^' k y,j 



+\^ k (-X\, - X k f + X l hk ) + *jM>j\ kl »'] (45) 



a 2 a 2 



IV'vI/ - 2 i'VM'p + X %i - 4x li *,ii - Ix'^l'., + -VM',^ . (46) 
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The evolution equations for the peculiar volume expansion scalar and for the shear im- 
mediately follow from Eqs. ©, © and (JJJ) 

l„o 3 o 1 



9 + 3H6 + ^6 2 + j-a 2 



4a 2 



ft, 



-i ( 7? 



-7^ 



(47) 



Replacing our expressions for 9, x l j, T^j an d 7£ in terms of the metric coefficients and 
retaining only terms containing up to four spatial derivatives, one obtains differential equa- 
tions for \I/ and Xiji namely 



kl , 



+ (x 4 %-ix fc W\)+3*X 



J a 2 



^, - -TV?, 



(48) 



(49) 



These equations can be solved iteratively. With two gradients only, the conformal Ricci 
tensor TV j and scalar 1Z coincide with their initial values, i.e., with the curvature of the 
seed conformal metric e~ 2 ^™5ij. Up to two gradients we obtain 



r is w v# . 



vV - ~ (v^) 

o 



- 



3 \a J \H 
where D i j = frdj - \V 2 5 l r 



,10v/3 



(50) 
(51) 



From Eq. (jBTIj) we obtain the volume expansion scalar and the shear tensor: 



1 fa\- 1 / 2 / 2 



He 



,10^/3 



3 Veto 

3 Va / V-Hq 



o 



(52) 
(53) 



Let us first explain how the renormalization group method works at the level of two 
gradients starting from the solution of Eq. (JHOj) written in the form 



e 2 * in e {a — a-, 



(54) 



18 KclqJ \Hq, 

where e«1. By taking the limit e -C 1 we can isolate the long-wavelength part of in 



other words, \EV 
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The constant ^ m represents the value of the gravitational potential at some initial instant 
of time when the scale factor is a- m . We regularize the 0(e) secular term by introducing an 
arbitrary "scale factor" \i and a renormalized constant ^ lVL = ^n{fi) + e <5\I/(/i, a in ). If we 
split the term (a — ai n ) into (a — fi + \i — ai n ), then to first order in e 

* = + e5^(/i, a in ) +ee 2 ^ R ^\a -// + //- a in ). (55) 

The counterterm is determined in such a way to absorb the (fi — ai n )-dependent term in 
the gravitational potential 

6V{n, a in ) + e 2 **^ - a in ) = 0. (56) 
This defines the renormalization-group transformation 

VrQjl) = *in + ee 2$aW (/i - a in ), (57) 
and the renormalization group equation 

^*(A0 =ee «Mri (58 ) 

The solution of Eq. (JH%|) is 

= -~ln(c 2 - 2 e/i), (59) 

where c 2 = e ~ 10 ^/ 3 is the constant of integration. Equating // to the generic scale factor we 
find that the renormalized improved solution for the gravitational potential at the level of 
two gradients is given by 



5 1 

* = * R Gu = a) = -<p - - In 



(60) 



9 Vao/ \HqJ V 6 

If expanded up to two gradients, this solution coincides with Eq. (|5L)p. Since by con- 
struction one should take the long-wavelength part of the argument of the logarithm in the 
solution of Eq. (|60p. it is easy to see that the latter matches Eq. (JSfiJl expanded up to two 
gradients. Indeed, write Eq. as 

*«<*■*> -*'<**J = -H££^)- (61) 

Inserting Eq. (J5Uj) into Eq. (pTTj) and expanding up to two gradients using the fact that 
^(x, ti n ) = 5(/?/3, we obtain 

1 (Id e- 3 ^(l-3e w ^ 3 ea] d 3 
^(x,t)-^( x ,t in ) — — - In ' ; 



3 \ J D e- 3 *in d 3 x 
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J D e- 3 *in d 3 x J 



6 



,(62) 

Any 

which coincides with Eq. ([60)1 . Notice, in particular, that Eq. (|6U|) differs from the toy 

gravitational potential adopted by Hirata and Seljak [26]. 

Let us compute the corresponding deceleration parameter 

H D H-iju 
V = —FT2- 1 = — / r-2, (63) 

where we have used the fact that H D = H — \£^. Inserting Eq. (|6l)|) into Eq. (|63|). we find 
that at late times 

q ~ \ ■ I ~ 1 = 0, (64) 
i.e., the deceleration parameter tends to zero. This result confirms our expectation that 
at the (resummed) lowest order in the gradient expansion, the Universe turns out to be 
curvature-dominated at late times, which is equivalent to a Universe with effective equation 
of state w = —1/3. What about the resummation of the long- wavelength perturbations at 
higher order in gradient terms? The curvature term is a series of gradients, and can be 
written in the form 

R=J2 e 2n ^c n a n -\ (65) 

n>l 

where c n — O (d 2n ) is a coefficient containing 2n gradients. Repeating the resummation 
procedure outlined for the case of two gradients, one can easily show that at any given order 
n the renormalized solution reads 

tffi(a) ~-^ln(l-2nc n a"). (66) 

Since at very late times ^ r ~ — |ln a, the corresponding deceleration parameter at late 
times goes like in Eq. (|64|). We conclude that at any order in the gradients, at late times 
the effect of the resummation of the long-wavelength perturbations is simply to generate 
a curvature term. This conclusion may be obtained also by inspecting Eq. (f63j) after the 
constant of integration \l/i n has been promoted to the renormalized quantity ^r. Each term 
in the series gives a contribution to R which scales as a~ 3 ~ t~ 2 . If only long- wavelength 
perturbations were present, the true scale factor a D would scale like 

a D = ae^ lna = a^ 2 , (67) 
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and (R) d would scale like a^> 2 . Therefore, if only long- wavelength perturbations were present, 
at large times and at any order in the gradients the line-element would take the form of a 
curvature dominated Universe, with hij ~ t 2 CV,(x), where CV,(x) is a function of spatial 
coordinates only. 

In summary, super-Hubble perturbations cannot be distinguished from the background 
for local observers. Thus a Universe which is pure matter and has only super-Hubble per- 
turbations, looks like a FRW universe to the local observer. Even if we started with a flat 
Universe plus perturbations, it is clear that the local observer will interpret what she/he 
sees as a FRW model with curvature (it would need a fine tuning to have k=0 within the 
Hubble patch). Now, as there is only matter and curvature in that model, the curvature 
will eventually dominate at late times, as a (open) non-flat matter Universe is dominated 
by curvature at late times. 

B. The effect of sub-Hubble modes 

Dealing with the backreaction of sub-Hubble perturbations, and therefore attacking the 
issue of the cross-talk between the sub-Hubble modes and the homogeneous mode, is more 
difficult than dealing with the super-Hubble modes because, as we shall see, the gradient 
expansion displays an instability of the perturbative series. 

In the effective Friedmann description of the inhomogeneous Universe one wishes to com- 
pute the typical value of the local observables averaged over the comoving volume D. By 
that we mean the ensemble average of such a volume average. The cosmological perturba- 
tions are treated as variables that take random values over different realizations of volumes 
D. In other words, we calculate the typical value of a quantity for a region of given size 
as the statistical mean over many different similar regions. This typical value is generically 
accompanied by a variance. If the size of the comoving volume D is much smaller than the 
global inflationary volume, then we can imagine placing this volume in random locations 
within a region whose size is much bigger than the size of D. By the ergodic property, 
this is equivalent to taking random samples of the ensemble for a fixed location of the box. 
In other words, one can replace the expected value of a given quantity averaged over a 
given comoving volume D with the ensemble average of the volume average, denoted by 
(■ • •) D . Since we are interested in the role of sub-Hubble perturbations which cause a tiny 
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variation of the value of the gravitational potential from one Hubble patch to another, the 
variance of the local mean observables is small. Under these circumstances, we can safely 
replace the spatial average with the ensemble average. This automatically implies that the 
perturbations which contribute to the effective dynamics are no longer restricted to receive 
contributions peaked at modes comparable to the Hubble-size (technically, this means that 
the window filter function defining the size of the comoving volume D plays no role) and 
therefore can be much bigger than of order 10~ 5 (or powers of it) 6 . 

Let us start by considering the lowest order in a gradient expansion, i.e., keeping only 
two spatial derivatives. The mean local curvature will be non- vanishing, but Qd will be zero 
at this order, as \I/ contains at least two spatial derivatives. In such a case, the integrability 
relation Eq. (J25|) immediately shows that the only consistent solution is 

(R) D oc a~ D 2 , (68) 

i.e., the effect of sub-Hubble perturbations at this order is to generate a standard curvature- 
like term in the effective Friedmann equations, scaling as the inverse square of the scale 
factor. This simple result holds at any order in perturbation theory (provided that one keeps 
only two spatial derivatives) and represents a straightforward extension of what found in 
Ref. [44] where it was shown that to second order in spatial gradients and in the gravitational 
potential, cosmological perturbations amount only to a renormalization of the local spatial 
curvature (this result valid up to two gradients was though improperly applied to the findings 
of Ref. P| , where more than two spatial derivatives were included, for instance through the 
physical redshift; this point was also noticed in Ref. y|). The result of Eq. (J58J) is reminiscent 
of the so-called vacuole model (see, e.g., Ref. |3|)- 7 Consider indeed a spherical region of a 
perfectly uniform Universe. Suppose that the matter inside that spherical region is squeezed 
into a smaller uniform spherical distribution with higher density. By mass conservation there 
will be a region in between the overdense sphere and the external Universe that is completely 
empty. Einstein's equation are exactly solvable for this situation in terms of the Tolman- 

6 It is important to stress that although the evolution of the kinematical backreaction and the mean curva- 
ture are obtained for averaged fields restricted to the domain D, the solutions to the averaged equations 
are actually influenced by inhomogeneities outside the domain D too, since the initial data are to be 
constructed non- locally and so take the fields on the whole Cauchy hypersurface into account. We thank 
T. Buchert for discussions on this issue. 

7 We thank S. Carroll for correspondence on this issue. 
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Bondi metric. The outside FRW Universe is totally unaffected by such a rearrangement. By 
Birkhoff 's theorem, the empty shell will be described by the Schwarzschild metric. Finally, 
the interior region will behave like a homogeneous and isotropic FRW Universe, but with 
different values of the cosmological parameters. These parameters will exactly obey the 
conventional Friedmann equations, and someone who lived inside there would have no way 
of telling that those parameters did not describe the entire Universe. 

However, the situation changes if we consistently account for higher-order derivative 
terms both in (R)d and in Q D . The lowest non-zero contribution to contains four 
spatial gradients and goes like a 2 H 2 oc a -1 ; this corresponds to a similar term with four 
gradients in the mean spatial curvature. 

Let us further elaborate on these findings. In all generality one can write 

oo 

Qd = Y,Qna n - 3 

n=2 

OO 

(R)d = E^« n " 3 > ( 69 ) 

n=l 

where q n and r n are expansion coefficients containing 2n spatial gradients. (Note that qi = 0, 
as starts from 4 gradients.) One may wonder about the actual range of validity of the 
gradient-expansion technique. At first sight it might appear to be valid only to describe 
inhomogeneities on super-Hubble scales, i.e., for comoving wave-numbers k ^ aH. However, 
this is not really the case! As one can easily check, terms of order n in the expansion, i.e., 
terms with 2n gradients, contain the peculiar gravitational potential (p to power m with 
2n > m > n. The dominant contribution at each order n {i.e., with 2n gradients) is 
Newtonian, i.e., coming from terms of the type (d 2 <p) n . However, these terms both in Qd 
and (R) d sum up to produce negligible surface terms when averaged over a large volume, 
so that the leading terms become the first post-Newtonian ones, i.e., those proportional to 
(d 2 (p) n ~ l (d(p) 2 . In other words the expansion is shielded from the effect of the Newtonian 
terms, which could in principle be almost arbitrarily large, by the volume averaging. The 
same protection mechanism, however, does not apply to the non-Newtonian terms in the 
expansion, simply because they cannot be recast as surface terms. This simple reasoning 
immediately leads to the conclusion that the actual limit of validity of our expansion at order 
n, is set by (k/aH) 2n (p n+1 < 1. Because of the nearly-Gaussian nature of our inflationary 
seed ip, it is clear that the lowest-order term able to produce a big contribution to Qd and 
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(R) d appears for n = 3, i.e., a term with six gradients. The importance of the six-derivative 



post-Newtonian terms has indeed been stressed also by Notari 
fourth-order moment of ip of the type 



It is a disconnected 



(W) 2 /^o 4 )((V^) 2 /^o 2 >, (70) 

having assumed that the spatial average coincides with the ensemble average. At this level 
an instability of the perturbative expansion is produced by the combination of the small post- 
Newtonian term ((V \p) 2 / Hq) (of order 10 -5 ) with the Newtonian term ((V 2 <^) 2 / 'Hq) , which 
can be almost arbitrarily large [6(, due to the logarithmic dependence on the ultraviolet 
cut-off (for a scale-invariant spectrum and cold dark matter transfer function). 8 

It is important to stress that the six derivative terms give a contribution to Qd which 
scales like a?H 2 = constant; similarly the six-gradients contribution to the smoothed curva- 
ture scales like a 2 /a 2 = constant. So these terms give rise to a sizeable effective cosmological 
constant-like term in our local Friedmann equations. In order to estimate correctly the 
six-gradients terms one needs the metric coefficients \I> and \ l j U P to four gradients (whose 
explicit expressions are given in the Appendix). 

The existence of a large contribution at six gradients, however, suggests that higher-order 
gradient terms will similarly lead to large corrections to the FRW background expansion rate. 
This is indeed the case. In the large-n limit there will be large contributions coming from 
perturbations in the quasi-linear regime (|5fw| ~ !)• These generic conclusions, however, 
also tell us that stopping the expansion at six gradients would be completely arbitrary and 
that, in any case, the perturbative approach cries for a more refined treatment than simply 
counting powers of the scale factor as done in Ref . |l2j . The existence of large corrections to 
the background should be taken strictly as evidence for an instability of the FRW background 
caused by nonlinear structure formation in the Universe. The actual quantitative evaluation 
of their effect on the expansion rate of the Universe would however require a truly non- 
perturbative approach, which is clearly beyond the aim of this paper. 

Connected to this fact is a technical obstacle in extending the validity of the gradient 
expansion to late times and/or to the nonlinear regime. This comes from the fact that 



Unlike the standard perturbative approach, we need not require small matter density fluctuations 
\7 2 (P/Hq <C 1, in our approach. This is because in the evaluation of mean obscrvables, powers of V 2 ^/^ 
either give rise to tiny surface terms or get multiplied by the small post-Newtonian term (W(p) 2 /Hq. 
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the metric determinant may become negative, indicating an internal inconsistency of the 



approximation. In Ref. [42j the problem is solved by using an "improved" approximation 
scheme which expresses the metric as a "square;" this choice guarantees non-negativity and 
leads to a GR extension of the classical Zel'dovich approximation of Newtonian theory. It 
is then shown that with suitable choice of the initial seed, such an improved approximation 
provides an excellent match to an exact inhomogeneous solution of Einstein's field equations, 
the so called Szekeres metric which describes locally axisymmetric (pancake) collapse 
of irrotational dust. Alternatively, exploiting the non-perturbative continuity equation, (1 + 
^frw) 1 — Ji* dtO, one can easily convince her/himself that the determinant of the metric 
is always well-defined. 

In order to take one step forward in the gradient-expansion approach, we will use the 
same renormalization group technique previously applied to deal with the backreaction of 
super-Hubble modes. Let us start by dealing with the case of two gradients. Can we apply 
the renormalization technique to the case of sub-Hubble perturbations up to two gradients? 
The answer is yes, since the spatial averages of objects like V 2 (p/Hq and (Vip) 2 /Hq can 
be replaced by the corresponding ensemble averages and therefore are tiny (of the order of 
1CT 5 ). The renormalized growing solutions at two gradients reads therefore 9 



5 1 , 

* = -ip In 

3^ 2 



(71) 



The next step consists in solving for the cosmological perturbations at four gradients. The 
equations of motion for the gravitational potential \l/ and for x l j at four gradients are given 
by Eqs. ()48|) and where the sources in the right-hand-side are computed inserting the 
solutions of Eq. (f7Tj) . 

Upon defining the coefficient 

^ = I e 10 ^ 3 fvV-|(V^) 2 V (72) 



9 V 6 



9 The long-wavelength part of the factor e 5v ( x ^ 3 , associated with each spatial gradient, can be re-absorbed 
by a redefinition of the spatial coordinates, as noticed in Ref. |2fij . and does not play any role when eval- 
uating ensemble averages as well as in the backreaction problem (we thank M. Porrati for correspondence 
on this issue); the small-wavelength part, on the other hand, can be expanded as (1 + 5ip s /3) ~ 1 because 
(p s ~ 10 -5 . Nonetheless, we prefer to show them explicitly because they provide the initial condition C- ln 
for the renormalization approach. 
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the matrix 

= (tfdjtp - ) , (73) 



j 

and the traceless matrix 

Ci _ ^ 10^/3 



6 



(74) 



the growing solution at four gradients for the gravitational potential assumes the form 



$ ~ — — Tr 
24 

+ 

36 S 



in 1 1 - 1 — (i) 2 s\ mfi-m 



ao/ \HqJ J \ \ao/ \Hq 



(75) 



A similar solution can be obtained for \ % • at four gradients. At late times the solution for 
the gravitational potential grows like (In a) 2 . A renormalization procedure can be applied 
to the solutions at four gradients because there are still "small" perturbative terms at hand, 
for instance terms like (Vy?) 4 whose spatial average is small. This amounts to saying that 
one has to take the solution of Eq. (|75jh expand the arguments of the logarithms and apply 
the renormalization procedure described previously at second order in the "perturbative 
parameter" £. The renormalized solution for the gravitational potential will grow like (In a) 
at large times. 

The lesson to learn from this computation is that, if we proceed further and go to six 
gradients, the unrenormalized gravitational potential, as well as x % ji wm grow like (In a) 3 . 
This is surely a step forward compared to the simple counting of powers of the scale factor 
which predicts that, at six gradients, the gravitational potential should grow like a 3 . How- 
ever, at this stage the renormalization procedure fails because it involves terms like the one 
in Eq. (|7U|) . which may be easily of order unity. Even the resummed perturbative expansion 
shows an instability produced by the combination of post-Newtonian and Newtonian terms; 
solutions with 2n gradients are expected to behave like (lna) n . If taken at face value, such 
a time-behavior of the gravitational potential would lead to acceleration of our local Hubble 
patch. To put this indication on firmer grounds (or to disprove it), however, one should go 
beyond the perturbative approach adopted in this paper. Our result may spur the efforts 
toward the search for a nonperturbative description of the dynamics of the system which 
would account for combinations of large Newtonian and small post-Newtonian terms. 

As a concluding remark of this subsection, we address a common objection to the use of 
the synchronous and comoving gauge in addressing the backreaction problem, namely that 
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the occurrence of shell-crossing singularities (caustics) in the evolution of collisionless fluids 
might prevent the analysis to be carried over into the fully non-linear regime. We would 
like to point out that the instability we find in the gradient expansion is unrelated to shell- 
crossing singularities. This can be immediately appreciated by noting that: i) shell-crossing 
instabilities imply the emergence of divergent gradients terms, while our instability shows up 
through an infinite number of finite gradient terms; ii) shell crossing is well known to lead 
to an infinite Newtonian term, while our effect involves a tiny Newtonian term. It should 
also be stressed that the occurrence of caustics does not represent a serious limitation of our 
approach; indeed, the very fact that caustics only carry a small amount of mass implies that 
they can be easily smeared over a finite region out in such a way that their presence does 
not affect the mean expansion rate of the Universe. 

For the sake of completeness, in the next subsection we will address the problem at hand 
within the commonly used weak- field approximation in the Poisson gauge. 



C. The backreaction in the weak-field approximation 



So far, in evaluating the effect of backreaction we have been making use of a perturbative 
approach in which non-linear dynamical quantities are explicitly expressed in terms of the 
inflationary seed perturbation (p and its spatial derivatives. This is the reason why higher 
and higher gradients of (p appear in our results. The same conclusion would hold also in 
different gauges as well as by using different perturbative schemes. 

One might however wonder whether the back-reaction problem can be approached directly 
in terms of non-linearly evolved variables. Related to this issue is the gauge choice. Non- 
perturbative approaches are indeed possible both in the comoving gauge adopted so far (see 



Ref. 



25|). and in the more commonly used Poisson gauge 47]. Working out the effects of 



back-reaction in a non-comoving gauge is indeed fully legitimate, provided a well-defined 
space-time splitting is performed, e.g., by means of the ADM approach 2^. We will here 
only sketch how back-reaction effects can be evaluated in the Poisson gauge, leaving to a 
subsequent paper a more detailed and quantitative analysis of the problem. 
The line-element of the Poisson gauge reads (see, e.g., Ref. 3|) 



ds 2 



! (r) {- (1 + 2<p P ) dr 2 - 2V i drdx i + [(1 - 2^ P ) <% + h$p] fo^dx?} 



(76) 
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where r is the conformal time and a(r) oc r 2 is the FRW background scale-factor for our 
irrotational dust source. It is important to stress that this line-element is meant to include 
perturbative terms of any order around the FRW background. The quantities Vi are pure 
vectors, i.e., they are divergenceless, d % Vi = 0, while h\j represent traceless and transverse 
(i.e., pure tensor) modes, h^ 1 = d l h^ = (spatial indices are raised by the Kronecker 
symbol). Vector and tensor metric modes are, respectively, of 0(l/c 3 ) and (9(l/c 4 ). To 
leading order in powers of 1/c the above line-element is known to take the well-known 
weak- field form 

ds 2 = a 2 (r) [-(1 + 2(f) P )dT 2 + (1 - 2i)p) d^d^dx^ , (77) 

where the scalars <pp and ipp are both 0(1/ c 2 ) and 4>p = ipp = $at/c 2 ; the Newtonian 
gravitational potential $at is related to density fluctuations 5p by the cosmological Poisson 
equation V 2 $at = 4irGa 2 5p. It is easy to realize that this form is accurate enough to describe 
structure formation within the Hubble radius as long as the considered wavelengths are much 



larger than the Schwarzschild radius of collapsing bodies 49]. 

The crucial point is that the kinematical back-reaction will contain the relevant term 
22| (N 2 Q 2 )d, where = (u^ being the fluid four- velocity). Here N is the inhomoge- 



neous lapse function needed to express the Poisson-gauge coordinate time tp = J dr a(r) 
as a function of the proper time t of comoving observers. This issue was already pointed 
out in Ref. |Q| where an approximate explicit expression for N was given (a second-order 
perturbative expression can be found in Ref. 0|); a term of the type (V$at) 2 appears ex- 



plicitly in N. What is important for us here is that Qp> will clearly display the same type of 
post-Newtonian (hence non-total derivative) terms which were found in the comoving gauge 
using our gradient expansion, namely terms of the type 



V 2 $„j (V$ w ) 2 ^, (78) 

where <& v is the velocity potential, which coincides (up to a sign) with the gravitational 
potential on linear scales; more generality <$> v and $at are connected by a cosmological 



Bernoulli equation (see, e.g., Ref. 25]). Similar terms appear in the mean curvature when 
projecting onto the comoving observer frame. We stress again that the terms of the type (J78|) 
appear only when considering the correct effective description of the average dynamics which 
has to include the kinematical backreaction term. Notice that this does not amount to saying 
that post-Newtonian effects are relevant in the dynamical evolution of the gravitational 
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and velocity potentials themselves. Indeed the expression (J75J) requires evaluation of the 
generally non-linear potentials and $at which may be readily obtained through the use 
of standard iV-body simulations. Owing to the non-linear (hence non-Gaussian) nature of 
the potentials and $tv ; the average (J78|) contains both a disconnected term, as in our 
previous treatment, and a non-zero connected four-point moment which is dominated by 
mildly non-linear scales, of order a few Mpc. 

Contrary to what happens in the synchronous gauge when the result is expressed in terms 
of the initial seed tp, in the weak-field approximation the number of gradients is expected 
to be finite and the complexity of the problem resides in the non-perturbative evaluation of 
the evolved potentials $„ and $at. 

It is interesting to note that the combination in Eq. (J78|) provides a contribution to Qe> 
which is of the order of H 2 and, using the linear dependence, nearly constant in time. 



IV. CONCLUSIONS 



The most astonishing recent observational result in cosmology is the indication that our 
Universe is presently undergoing a phase of accelerated expansion. One possible explanation 
of the observations is that the Universe is homogeneously filled with a fluid with negative 
pressure that counteracts the attractive gravitational force of matter fields. Another possible 
explanation is a modification of GR on large distance scales. 



In this paper we have elaborated on the alternative idea that the backreaction of cos- 
logical pe—ons may caase the co smi c acceleration Q Q fl Q fl fl 
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1 Q. 



Following Buchert 21, 22], we have provided the effective Friedmann equations describing 
an inhomogeneous Universe after smoothing. The effective dynamics is governed by two 
scalars: the so-called kinematical backreaction Qo and the mean spatial curvature (R)d- 
They enter in the expression for the effective energy density and pressure in the Friedmann 
equations governing the mean evolution of a local domain D. For positive Qd, accelera- 
tion in our local Hubble patch may be attained despite the fact that fluid elements cannot 
individually undergo accelerated expansion. Indeed, the very fact that the smoothing pro- 
cess does not commute with the time evolution invalidates the no-go theorem, which states 
that there can be no acceleration in our local Hubble patch if the Universe only contains 
irrotational dust. 
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Through the renormalization group technique, we have then shown that super-Hubble 
modes can be resummed at any order in perturbation theory yielding a local curvature 
term ~ a^ 2 at large times. We then turned our attention to the backreaction originating 
from modes within our Hubble radius, studying perturbatively their time-behavior. In this 
case our findings indicate that an instability occurs in the perturbative expansion, which 
may be not taken care of by the renormalization group procedure since terms of the form 
H 2 (5p RW (v/c) 2 ) (where v is the peculiar velocity) start appearing both in Q D and in the 
mean spatial curvature. Such terms are not as small as order 10 _5 if 2 ; on the contrary 
the averaging procedure allows the combination of post-Newtonian and Newtonian terms 
to acquire values of order H 2 . Since the perturbation approach breaks down, we may not 
predict on firm grounds that backreaction is responsible for the present-day acceleration of 
the Universe. However, it is intriguing that such an instability shows up only recently in the 
evolution of the Universe and that this picture is further supported by a very general result; 
as shown explicitely by Buchert et al. 24|, even a tiny back-reaction term can drive the 
cosmological parameters on the averaging domain far away from their global values of the 
standard FRW model, thus modifying the global expansion history of the Universe. Other 
aspects of the scenario discussed in this paper, such as the dynamics of perturbations on 
observable scales, will be the subject of a forthcoming publication. 



APPENDIX: FOURTH-ORDER GRADIENT-EXPANSION APPROXIMATION 
TO THE SOLUTION OF EINSTEIN'S FIELD EQUATIONS 



For completeness, we give here the explicit expression for \1/ and x% up to four gradients 



fwe refer the reader to Ref. 



for the detailed derivation of these results). We have 
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+ ™(Vv,) 2 (^V.-^(V^)}. (A.2) 

One can verify that these solutions satisfy the energy constraint and the momentum 
constraint up to the relevant number of gradients. It is also important to stress that these 
expressions reproduce the perturbative second-order metric (see, e.g., Ref. p) when only 
terms up to second order in ip are kept. 
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